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Abstract. A hypercomputation model named Infinite Square Well Hy- 
percomputation Model (ISWHM) is built from quantum computation. 
This model is inspired by the model proposed by Tien D. Kieu \£\ and 
solves an Turing-incomputable problem. For the proposed model and 
problem, a simulation of its behavior is made. Furthermore, it is demon- 
strated that ISWHM is a universal quantum computation model. 

1 Introduction 

In general terms, a quantum computation model is built based on a physical refe- 
rent which determines the mathematical and computational characteristics of the 
quantum algorithms built on such model. These algorithms, in turn, determine 
the kinds of problems that can be solved under the quantum computation model. 

The first quantum computation model considered as a hypercomputation 
model was the one proposed by Tien D. Kieu [J having the following characteris- 
tics: (i) Physical referent: quantum harmonic oscillator, (ii) Turing-incomputable 
problem that solves: Hubert's tenth problem. 

This article presents the construction of a hypercomputation model named 
Infinite Square Well Hypercomputation Model (ISWHM). This model is based on 
the one proposed by Kieu, but we have selected the infinite square well instead 
of the quantum harmonic oscillator as the physical referent. This change conveys 
a change in the Turing-incomputable problem to solve, and instead of solving 
Hilbert's tenth problem it solves a problem equivalent to it. A simulation of 
its behavior is carried out for the proposed model and problem. Besides, it is 
demonstrated that ISWHM is a universal quantum computation model. 

2 Hypercomputation 

While the idea of an absolute computability, detached from logical, mathematical, 
physical or biological theories is hard to support nowadays, the idea of a relative 
computability has progressively gained supporters, as shown by the establishment 
of an academic community around the topic (www.hypercomputation.net). 

The term 'hypercomputer' or 'hypermachine' denotes any data processing 
device (theoretical, potentially realizable or that can be implemented) capable 



of carrying out tasks that cannot be performed by a Turing machine A 
hypermachinc, also capable of simulating a universal Turing machine is called a 
super-(Turing machine), otherwise it is called a non-(Turing machine) 

At first glance, we could think that the possibility of a hypercomputation 
model would be a refutation of the widely accepted Church- Turing thesis, which 
identifies the naturally calculable functions with the Turing-computable func- 
tions. But actually the existence of hypercomputation models refute the M 
thesis, which identifies the functions calculable by a machine with the Turing- 
computable functions 0]. On the other hand, notwithstanding the proliferation 
of theoretical hypercomputation models [S] , the possibility of the real construc- 
tion of a hypermachine keeps being controversial and under analysis. 

3 The Turing-Incomputable Problem 

The demonstration that ISWHM is a hypercomputation model will be done by 
demonstrating that it can solve an Turing-incomputable problem. 

Let IN be the set of non-negative integers, 7Z the set of the whole numbers, 
7Z + the set of the positive integers, and 7Z +2 the set of the perfect square 
numbers. A Diophantine equation is one of the form 



where D is a polynomial with integer coefficients. In present terminology, Hilbert's 
tenth problem may be paraphrased as: Given a Diophantine equation of type 
we should build a procedure to determine whether or not this equation has a so- 
lution in ZZ. From the concluding results gotten by Matiyasevich, we know that, 
in the general case, this problem is algorithmically insolvable or more precisely, 
it is Turing-incomputable. 

Having into account that any non-negative integer number may be expressed 
as the addition of the squares of four integer numbers (the four squares theo- 
rem, due to Lagrange), and that any system of a finite number of Diophantine 
equations has a solution in ZZ^ if and only if a certain Diophantine equation 
associated to the system has a solution in 2Z\ the Hilbert's tenth problem may 
be reduced to search a solution in IN 0. 

Due to the Hamiltonian operator J5J used in the description of ISWHM, 
the decision problem associated with this model of computation consists in de- 
termining whether or not a Diophantine equation kind JQ) has a solution in 
-2T +2 . However, this problem is equivalent to determining whether or not the 
Diophantine equations system 



has a solution in 2Z, and this problem is equivalent to Hilbert's tenth problem. 
Therefore, the decision problem associated with ISWHM is a Turing-incomputable 
one. 
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4 The Physical Referent 



The energy levels E n , for a particle with mass m in a infinite square well with 
length L, are given by [7] 
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These energy levels are the eigenvalues of the Hamiltonian operator of the 
system denoted by H, for which correspond as eigenvectors the states denoted 
by the kets | n), where n <G 7Z^ . The Hilbert space for the infinite square well, 
denoted 1H C , is infinite dimensional, and its orthonormal and complete canonical 
base is the countably infinite set {| n )}„ 6 ^+- I n such base, H is diagonal and 
has the following form 



H\n) = ^\n). (2) 



If matrix M is defined as 



then @ may be written 
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M = — = H , 



M\n)=nf\n) . (3) 

Equation © implies that matrix M has a diagonal form in its canonical 
base {\ n )} ne zZ + an< ^ ^ s com P onen ts are the 2Z +2 numbers. This canonical 
base satisfies the following completeness and orthonormality conditions 

oo 

^ | n) (n 1 = 1, (n \ m) = S n , m , 
n=0 

where 11 is the identity matrix in Hilbert space IH C . 

With this mathematical description of the infinite square well we may proceed 
to the extension of Kieu's model from the quantum harmonic oscillator towards 
the infinite square well thus to obtain the ISWHM. 



5 The Algorithm 

If we have a Diophantine equation of type Q , the strategy to determine whether 
or not it has a solution in 2Z +2 , inspired in the strategy followed by Kieu, consists 
in translating equation D in k variables into a codifying Hamiltonian operator 
Hd represented by a countably infinite dimensional matrix which acts in Hilbert 
space Hf k = H c ® • • • ® H c , and then restate the problem in terms of the Ho 

k— times 

properties pQ. 



Particularly, each of the k variables is substituted in D by the operator M 
given by © and we obtain the Hamiltonian Ho given by 

H D = D(M 1 ,M 2 ,...,M k ) 2 , 

and therefore, equation D has at least a solution in 2Z +2 , if and only if the 
eigenvalue Eq associated to the Ho fundamental state is zero. 

The problem of determining the Eq eigenvalue may be solved by the appli- 
cation of the adiabatic theorem, as Kieu does in the harmonic oscillator case. 
For this doing a universal interaction Hamiltonian operator Hi (valid for all 
Diophantine equations) is introduced, which describes an operator that acts in 
Hilbert space H® fc and has a eigenvector known in such a space that corresponds 
to the zero eigenvalue. 

The Hamiltonian operator H^(t/T) that results from the convex overlapping 
of Ho and Hi, with adiabatic parameter t/T has the form 



H A (t/T) = (1 - t/T)Hi + (t/T)H D , 

where T is the total time of the adiabatic evolution and t S [0, T]. The Hamil- 
tonian H^{t/T) determines the solution of the Schrodinger discretized equation 
given by [8] 

where St is the time step. 

For a Diophantine equation D(x) — 0, based on Q the algorithm for ISWHM 
to determine whether or not it has a solution in ZZ +2 , has as inputs and outputs 
the ones indicated in Fig. ^ where for ip(t) = c n (t) \ n): 

1. P n {t) — |c„(t)| 2 represents the probabilities of the quantum states | n) for a 
time t. 

2. < n 2 (t) >= J2 n Pn(t)n 2 represents the expected value of the square random 
variable n that codifies the value of the unknown x in the equation D{x) = 0. 

3. Eo(t) represents the spectral flow of the eigenvalue associated to the funda- 
mental state of the Hamiltonian H A (t/T). 

According to the characteristics of the Hamiltonian H^(t/T), only a finite 
adiabatic evolution T time is necessary to obtain a state | no) which has a high 
probability to be the fundamental state of Ho |SJ. The | no) state is observed 
from the emergence of a dominant probability in the set of probabilities P n (T). 
Once the | no) state is obtained, its eigenvalue Eq{T) is determined. If Eq(T) = 0, 
the Diophantine equation D{x) = has a solution in TZ^ 2 and it is given by the 
value < n 2 (T) >. On the other hand, if E (T) ^ 0, then the equation D(x) = 
doesn't have a solution in ZZ, +2 . The above description should be generalized for 
a Diophantine equation of k variables. 
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Fig. 1. Algorithm for the ISWHM 



6 The Simulation 

The idea is to simulate the algorithm described in Fig. in such a way that the 
realization of the adiabatic theorem for a big but finite T time is observed. Since 
the algorithm operates upon a infinite dimensional Hilbert space, it is necessary 
to operate upon a truncated finite dimensional space to obtain the numeric 
solutions indicated by (J3J. The truncation level P means that we operate upon 
a space which canonical base is given by {| f ),..., | P)}- 

For the carried out simulations, ^(O) is assumed as uniform and normalized, 
Hi is obtained from the matrix of ones, and 5t = 1. Concretely, we carried out 
the simulation for the following simple Diophantine equations 



x- 16 = , 
x- 7 = , 
(x + l){y + 2) = 12 



(5) 
(6) 
(7) 



where (0 has a unique solution in ZZ, +2 , © doesn't have a solution in TZ^" 1 and 
(J2J has a unique solution in £Z +2 x 2Z +2 . 

For (JjJJ, with T w 2000 and P — 6, it was determined that it had a solution. 
Figure 2(a) points that rapidly (t w 15) the state | no) = | 4) is the candidate 



to be the fundamental state Hd. Figure 2(b) indicates that this state effectively 
has the energy eigenvalue E = 0, and Fig. 2(c) shows that the solution is given 
by x — 16. 

©, Fig. 3(a) indicates that the probability that the state | 3) be the 



For 



fundamental state of Ho can be approximated to 1 as much as we want. However, 
Fig. |3(b)| points that Eq ^ 0, therefore, © doesn't have a solution in 2Z +2 . 

For J7J, Fig. 4(a) shows that the state |no,mo) = 1 1 , 2} tends to be the 
fundamental state, when T — 1000 and P = 2, and Fig. 4(b) shows the contour 



curves for the codifying variables n and m at three different time moments. 
Figure 4(c) indicates that the equation solution is given by x = 1 and y = 4, 
supported by the fact that the spectral flow of eigenvalue Eo(t) is similar to the 
one illustrated in Fig. 2(b) 



7 The Universality 



The universality of a quantum computation model, that is, its capability to 
carry out any operation that is carried out by a Turing machine, is established 




Fig. 2. Simulation results for© with T = 2000 and P = 6 




Fig. 3. Simulation results for JBJ with T = 2000 and P = 6 




Fig. 4. Simulation results forJTJ with T = 1000 and P = 2 



by its capability to generate a set of quantum gates, such that any unitary 
transformation U(2 n ), that is, any quantum gate that operates upon n-qubits, 
can be approximated with sufficient exactness by a quantum circuit that is only 
made of a finite number of gates of this set. 

Henceforth, we will use the convention of superindexes over the operators 
(quantum gates) and over the states (qubits). This superindex will denote the 
dimension of Hilbert space upon which the operators act or upon which the 
states are defined. 

To demonstrate the universality of the ISWHM it is necessary and sufficient 
to construct three quantum gates: the controlled-NOT gate C^ OT , the phase 
gate fi^M and Hadamard gate H^ d |9I1U| . For gates C^ OT and ^/ 4 > its matrix 
representation and the transformation carried out upon the base states of Hilbert 
space on which they act, are given by 
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|0) 2 |0) 2 , 



[0 e ifJ > |l)2^ e W4|l)2 . 

In systems whose Hamiltonian operators are independent from the time, as 
it is the case of the infinite square well, the states of the system evolve according 
to the solution of Schrodinger equation of stationary states [7] 

\m) = u(t)\m) = e-* m \m) > m 

where H is the Hamiltonian of the system defined in J3J), U(t) is the evolution 
unitary operator and | ^(0)) = \n)°° is considered for the eigenstates of the 
infinite square well Hamiltonian. The U (t) matrix elements are given by 

U np (t) = exp f-i 2^ 2^ ) Snp ' ^ 

where 6 np is Kronecker delta. 

From different choices of the t parameter in (|llfl . and different qubits coding, 
it is possible to find the evolutions corresponding to each of the mentioned 
quantum gates. 



In order to build the phase gate ^^y 4 in the ISWHM, based on the normalized 
eigenvectors of the phase gate ^ 2 / 4 ©, we carry out a coding of the canonical 
base {| 0) 2 , | l) 2 } for a 1-qubit in the canonical base {| n )°°} ne 2Z + ' Si yen by 

I „, 2 codification 1lC o 

I o) ► 1 1) , 

I -^2 codification | ^ oo 

According to and for the time t — 2 ^ 4>, where 4> S [—if, 7r), we obtain 
a quantum gate r°°(<fi) whose components are given by 

r~(0) = [exp (-i<j>) exp (-z0(n 2 - 1))] S np , (12) 

where term exp (—i(f>) hasn't physical meaning. 

With H12|) it is possible to control, at least theoretically, the ISWHM, until 
the implementation of the transformation that accounts for the phase gate 
given by 



#~ /4 = r~(-7r/12) |iy» 

/ 9tT7 r, 2 / 7T \ \ 

U 
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On the other hand, in order to build the C^q T gate in the ISWHM, based 
on the normalized eigenvectors of the C^qt S a * e ©j we c °de the base 



{l oo,Moir,^-^,^^, 

for a 2-qubit in the canonical base {| n }°°} nG ^+ by means of 

I qq^4 codification^ | 2^°° 
,„ ,4 codification ,,oo 

I 01) ► I 4) , 

| 10) 4 + | ll) 4 codification | 6)°° + | 1)°° 
|11) 4 - |10) 4 codification | 6) °° - |1)°° 

According to (|ll|l . in t = we obtain the C^ OT gate given by 
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In a similar fashion, the if^ gate is built in the ISWHM by making a coding 
of the {| 0) 2 , 1 l) 2 } base in the {| n )° C } ne 2Z+ base, from the normalized eigen- 
vectors of the H% d Hadamard gate, and by selecting an adequate value of the t 
parameter in 

Although from a theoretical point of view we got the universality of the 
ISWHM, this model doesn't turn out to be adequate from its physical imple- 
mentation point of view due to the realization of the mentioned quantum gates 
in different times and due to the change in the coding of the base of the qubits 
for each gate. 



8 Conclusions 

It is quite surprising that a quantum computation model over such a simple phys- 
ical referent like the infinite square well have hypercomputation characteristics 
as it was demonstrated for the ISWHM. In fact, and according to the obtained 
universality results, the ISWHM constitutes itself (at least from a theoretical 
perspective) as a super- Turing computational model. 

On the other hand, the success obtained by choosing a physical referent 
different from the one selected by Kieu opens the possibility of obtaining new 
hypercomputation models supported on the quantum computation. Due to the 
nature of their theoretical construction, these models may be the best candidates 
to a possible physical realization or implementation of a hypercomputer. Be 
this possible, such hypercomputer would broaden the conception of the term 
'computable', not exclusively from the theoretical perspective but also from the 
practical perspective. 
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